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Abstract 

Several aspects of complex- valued potentials generating a real and positive spec- 
trum are discussed. In particular, we construct complex-valued potentials whose 
corresponding Schrodinger eigenvalue problem can be solved analytically. 

PACS: 03.65.-W 

1 Introduction 

Quantum systems characterized by non-hermitian Hamiltonians are of interest in several 
areas of theoretical physics. For example, in nuclear physics [0] one studies standard 
Schrodinger Hamiltonians with complex-valued potentials, which in this connection are 
called optical or average nuclear potentials. Non-hermitian interactions are also discussed 
in field theories, for example, when studying Lee- Yang zeros 0. Even in recent studies 
on localization-delocalization transitions in superconductors P] and in the theoretical 
description of defraction of atoms by standing light waves non-hermitian Hamiltonians 
are of interest. 

Recently, several authors [|, |^, |^, § have studied standard one-dimensional Schrodinger 
Hamiltonians with complex- valued potentials giving rise to a real energy spectrum. Some 
of this work ^ has been concentrated on the numerical study of parity and time-reversal 
(PT) invariant Hamiltonians, and it is believed that this invariance is a sufficient crite- 
rion for the reality and positivity of the spectrum. It is the aim of this note to construct 
via the Darboux method new complex potentials for which the corresponding eigenvalue 
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problem can be solved exactly. These potentials are not necessarily PT invariant but still 
give rise to a real and positive spectrum. 

In order to make this note self-contained and also for the purpose of setting up our 
notation we briefly review in the next section the basics of the Darboux method. Using 
this method we then derive a rather general class of complex potentials giving rise to 
the harmonic oscillator spectrum. In particular, we note that this class of potentials also 
contains members which are not PT invariant. In Section 4, we will discuss the Darboux 
approach for the class of potentials V{x) = —{ix)^/2, which has recently been discussed 
in some detail by Bender and Boettcher |^ |^. Finally, in Section 5 we discuss the exactly 
solvable potential V{x) = |exp{2ix}, which may be considered as a superposition of the 
Bender-Boettcher potentials. 



2 The Darboux Method 



In this section we briefly review the method of Darboux [|10], which relates the spectral 
properties of a pair of standard Schrodinger Hamiltonians 

i^l/2 = -^^ + Vl/2(x) (1) 

acting on the Hilbert space L^(M). Let us assume that for one of these Hamiltonians, 
say Hi, the spectral properties are exactly known. That is, its eigenvalues En and corre- 
sponding eigenfunctions (pn are known explicitly: 

HiMx) = EnMx) ■ (2) 

For simplicity we will assume that Hi has a purely discrete positive spectrum enumerated 
by n = 0, 1, 2, . . . such that < Eq < Ei < E2 < . . .. If we now postulate that there 
exists a linear operator A such that 

AHi = H2A (3) 
then the functions ipn '■= A(j)n 7^ are obviously eigenfunctions of the other Hamiltonian 

H2MX) = EnMx), (4) 

with the same eigenvalue En > 0. 

A rather general form for an intertwining operator obeying (|^) is given by 

fc=0 

where the : M — * C, k = 0,1, . . . , N — 1, are (at least twice) different iable functions to 
be determined via the condition (|^) and Jn is an arbitrary constant. 
The simplest non-trivial choice is iV = 

A = -^ + /(x), /:M^C. (6) 



^More general approaches for standard real- valued potentials can be found in |Tl| 
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Putting this as an ansatz into eq.(^ leads to (/' = d//dx etc.) 

(^2 -Vr + - [{V^ - V,)f + Vi - f"/2]l = . (7) 

As ^ and the unit operator 1 are hnearly independent, their coefficients have to vanish 
simultaneously. That is, we have the two conditions 

V2 = V,-f', {V2-V^)f + V;-f"/2 = 0. (8) 

The ffist equation expresses the potential V2 in terms of Vi and /, and thus allows to 
eliminate V2 from the second one leading, after an additional integration, to 

f' + f - 2Vi = -2e = const. (9) 

Substitution / = u'/u then results in the Schrodinger-like equation 

1 d^ 



2dx2 



+ Vi{x) ] u{x) = eu{x) (10) 



with in general complex integration constant e. We will, however, restrict ourselves to 
real e for reasons to be given below. Note that we do not require square-integrability of 
u. We also note that in terms of / the two potentials read 

V^/2{x) = lP{x)±y{x)+e, (11) 



that is, for real / the two Hamiltonians H1/2 are SUSY partners [1^ . 

The above approach may be used to construct also complex-valued potentials which 
generate a real- valued spectrum containing that of Hi. To be more explicit, we start 
with a given real potential Vi, solve eq. (|To|) and thus obtain a new (in general complex) 
potential V2, which is given by 

V2{x) = (^y -Vi{x) + 2e. (12) 



u{x) 

Note that even for real e the above potential may be complex by choosing a complex 
linear combination of the two fundamental solutions of (|iy). This linear combination is 
not arbitrary because u must not have zeros on the real line in order to lead to a well- 
defined Hamiltonian H2 on L^(]R). Note that for e > Eq any regular solution u G L^(]R) 
will have zeros. Hence, we will respect the condition e < Eq, which is still not sufficient 
for having a u without zeros. 

With the help of the intertwining relation (^ it is now easy to see that the strictly 
positive spectrum of Hi forms a subset of the complete spectrum of H2 with the corre- 
sponding eigenfunctions given by 

lljn{x) = CnA(t)n{x) = C„ ( '(p'^x) + 0„(x)^ , (13) 

V u[x) J 

where C„ stands for a normalization constant defined by |C„|~^ = {A(f)n\A<Pn) ■ Note that 
this normalization constant vanishes if ipn is not square integrable, implying that the 
corresponding eigenvalue does not belong to the spectrum of H2. This, however, happens 
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only for the square-integrable solutions of ([T0|), which we have eliminated via the condition 
e < Eq. 

Noting that from (g) follows Hi{A'<)* = {A^)*H2 where (A^)* = d/dx + /(x), we realize 
that the above set of eigenfunctions ( ]TB| ) is only complete on L^(]R)\ker (A^')*. Note that 
{A^y is a first-order differential operator and, therefore, the dimension of ker (A"'")* may 
not exceed unity. In other words, H2 may have one additional eigenvalue which is below 
Eq. Note that ker (A^)* is a one-dimensional subspace of L^(M) iff the differential equation 
[A'^Yips = has a solution in this Hilbert space. This solution can explicitly be given 

M^) = ^ (14) 
u{x) 

and it may easily be verified that it is an eigenfunction of H2 with eigenvalue e, which 
belongs to the spectrum of H2 iff ipe ^ i^^(ffi). 

As conclusion of this section we note that the Darboux method may be generalized such 
that one can construct complex potentials generating a real spectrum which is identical 
to that of a self-adjoint standard Schrodinger Hamiltonian Hi, special = {Eq,Ei, . . .}. 
In some cases, via an appropriate choice of the parameter e and the solution of ( ]TD| ) 
this complex potential may have an additional real eigenvalue e which lies below of the 
spectrum of Hi. In the next section we will demonstrate this approach for the case of the 
harmonic potential Vi{x) = \ x'^. 



3 Complex potentials generating a harmonic 
spectrum 

In this section we will consider as an explicit example the harmonic oscillator potential 
Vi{x) = (l/2)x^. The spectral properties of the corresponding Schrodinger Hamiltonian 
Hi are well known: 



Er, 



n + 



2 ' 



[X 



[v^2'^n!]-^/2^„(x) exp{-a;V2} 



(15) 



where Hn denotes the Hermite polynomial [|T^ of degree n G {0,1,2,...}. The most 



general solution of the Schrodinger-like equation (|iy) can be given in terms of confluent 

,x')+(3xiFii'^/l,x')] . 



hypergeometric functions |p!3 



uix) 



[a iFi 



l-2e 1 
4 ' 2' • 



(16) 



As u should not have any zero on the real line a must not vanish and, therefore, can be 
set equal to unity without loss of generality. Note that an overall factor in u is irrelevant 
for the relevant formulas (p!2D -(p!4D. From the general discussion of the last section we also 
have the condition e < Eq = 1/2. Finally, we have to determine for which values of the 



remaining parameter /? G C the general solution (|T6D has no zeros on the real line. These 
conditions can be extracted M, IT5[ from the asymptotic behaviour 



Uix] 



exp{x^/2} 



\x 



1/2+e 



Fi 



+ f3 



F 



■3-2e^ 



For P G 



47 ^ V 4 7 

a strictly positive solution is only possible if 



X 



\x\ 



+ 0{l/\x\) 



(17) 
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Violation of this condition results in singularities in V2. See ref. [[T^ , ^ for details. 
However, for f3 G C/M the solution u does not have any zero. That is, the allowed values 
for P are given by the complex plane with the exception of the two cuts ] — 00, — /3c(e)] and 
[j3c{s), oo[ on the negative and positive real line, respectively. In addition, from (^) we 
can read off that ip^ = C^/u is square integrable in the allowed ranges of the parameters. 
In other words, we have an additional eigenvalue e < 1/2 in 7^2- In Figure 1 we show the 
real part (a) as well as the imaginary part (b) of the potential V2 for —3 < e < 1/2 and 
P = i. We note here that for all /5 G iM the potential is invariant under PT transformations 
1^, that is, V2*(— x) = V2{x) and, therefore, also the Hamiltonian H2 has this symmetry. 
However, other complex values of (3, which are not purely imaginary, are also admissible 
and thus provide examples of non-PT invariant potentials generating the real harmonic 
oscillator spectrum with the additional ground-state energy e. 

In Figure 2 and 3 we give several examples for fixed e = —1/2 and typical values of 
j3. Note that for this particular value of e the solution ([16|) can be expressed in terms of 
the error function |1l3| 



uix] 



Erf (x) 



This example explicitly shows that for all 



/5,(-l/2) = 2/v/^^ 1.12^ 



/5eC\{]-oo,-^]U[^,oo[} 



(19) 



(20) 



the solution given in ([T9|) does not have zeros on the real axis, because its imaginary 
part vanishes only at a; = where its real part is obviously non-zero. Figure 2 shows the 
real (a) as well as the imaginary (b) part of V2 for Re /3 = 0.5 and Im /3 G [—2,2], thus 
crossing the real axis of the complex /3-plane in the allowed region. Figure 3 shows the 
same for Re /3 = 2 and Im /5 G [—1,0], that is, reaching the cut of the /?-plane, which is 
lying on the positive real axis, from below. Here a singular behaviour is clearly visible, 
when /3 approaches the cut. 

We conclude this section in summarizing the spectral properties of }i2- 



specif2 = {e-, |, §,§,...}, ^e(x) = Celu{x) 



exp{— x^/2} 



[v/?2"+%!(n + l/2-e)]i/2 



n+l 



x) + 



u\x) 



- X HJx) 



(21) 



where the normalization constants have been determined via0 = (A0„|y40„) = 

2{(j)n\Hi — = 2{En — e). We note again that the spectrum of H2 is real and bounded 
below for any finite value ofe < 1/2. This remains true even in the case of complex f] where 
H2 is neither self-adjoint nor PT invariant in general. Because of lim|^|^oo Ini V2{x) = 0, 
the eigenvalue problem of H2 is well-defined on L^(R). Let us also remark that follow- 
ing the approach of [Q we are able to define ladder operators for the non-hermitian 



Hamiltonian H2 closing a quadratic algebra. 



^Note that this is a priori only true for (3 E 
continue also to complex values. 



But as the result is independent of (3 we may analytically 
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4 Application to the Bender-Boettcher potentials 



The Darboux method reviewed in Section 2 is also apphcable to complex potentials. Here 
we will consider the class of potentials,^ N > 2, 

T4(x):=-l(ix)^, (22) 

which has recently been discussed extensively by Bender and Boettcher 0. These 
authors have shown (numerically) that for N > 2 the above potential generates a real 
and strictly positive spectrum on an appropriate Hilbert space, which is taken to be the 
linear vector space of square integrable functions on some contour in the lower complex 
x-plane. This contour has to lie in some sectors centered about anti-Stokes lines and 
bounded by Stokes lines for Re x ^ ±oo 0. For simplicity, we will assume that this 
contour, which is the one-dimensional configuration space of the problem, approaches the 
anti-Stokes hues, that is, 

X ( n 2n } 
lim ■; — 7 = exp < —1— ± 1 > . (23) 



For the present case the Schrodinger-like equation (|10|) cannot be solved analytically 
for arbitrary e. However, for the special case e = 0, which is below the ground-state 
energy of Hi, this equation reduces to Bessel's differential equation. In other words, for 
this special case we can express the solution of (|l^) in terms of modified Bessel functions 

m 

u{x) = [a I,{z) + p K,{z)] , (24) 

where we have set 

^:=^— (ix)(^+2)/2 (25) 

Note that the anti-Stokes lines in the x-plane are mapped onto the negative real line in 
the z-plane (this is the cut of the Bessel functions) such that 

Hm ^ = e±'^ . (26) 

Re x— utoo \z\ 

The new potential generated via the Darboux method now reads 

'alu-Az) - l3Ky_i{z) 



V,{x) = -\{\x)'^ 



2 



(27) 



ah{z)+ (3K^{z) 

and approaches for |x| — > oo asymptotically the original potential Vi given in (P2|), that 



is, V2{x) = Vi(a;)[l -|- 0(l/|a;|)]. As a consequence, the new potential V2 has the same 
Stokes and anti-Stokes lines as the original one and, therefore, both Hamiltonians can be 
defined on the same Hilbert space. In other words, the Darboux transformation leaves 
the Hilbert space invariant. 

We remark that for appropriate choices of the parameters a and (3 the Hamiltonian 
H2 has an additional eigenvalue e = 0, which is its ground-state energy. This, in fact, 
corresponds to a situation where SUSY is unbroken. Note that the Hamiltonians Hi 12 are 

■^Note that the special case N = 2 corresponds to the harmonic osciUator discussed in the previous 
section. 
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SUSY partners |]T2|. In the case of broken SUSY H2 has the same spectrum as Hi. By 
allowing for complex a and (3 it is also possible to construct potentials V2 which are not 
PT- invariant. Nevertheless, the spectrum is identically with that of Vi which, for N > 2 
is strictly positive and discrete Note that in the present context PT-invariance of a 
potential V means V{x) = V*{—x) for all x G M. 



5 Another exactly solvable complex potential 

In this section we will consider the eigenvalue problem associated with the complex po- 
tential V{x) = I exp{2ix}, which may be considered as a superposition of all Bender- 
Boettcher potentials. This potential is periodic, V{x + tt) = V{x), and, in some sense, 
simulates a proper regularized large N limit of (0). The corresponding eigenvalue prob- 
lem can be write in the form 



^"(x) + (2E-e2'^) = 



:28i 



which in turn is reducible to Bessel's differential equation . Pairwise linearly indepen- 
dent solutions are, for example, given by 



ipi^x) 



/2E 



J. 



H 



(2) 
/2E 



(e-) 



(e'^O for E > 
J-V2E (e''^') for E < 



(29) 



In the complex upper-half x-plane, that is Imx > 0, the normalizable solution for 
> is given by ip^. For E' < we have ips and ip^ as normalizable solutions in 
the sectors < argx < 7r/2 and 7r/2 < argx < vr, respectively. There are no further 
restrictions on the parameter E. Hence, the spectrum of the corresponding Hamiltonian 
is given by the complete real line and thus is not bounded from below. In other words, 
in the complex upper-half x-plane the potential under consideration does not model a 
physically relevant system. The same holds for x G M, where no normalizable solutions 
can be found. 

In the complex lower-half x-plane, Imx < 0, the argument of the Bessel functions 
z := e^^ = el^™^le'^®^ becomes infinitely large for Imx —00. In order to find the 
physically acceptable sectors in the lower-half plane we consider the asymptotic behaviour 
of the pair ipi and ip2, 



iji{x) 



-e-'^/^exp 

TT 



exp 



2 A/i 
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+ 0(1^)), 

(l + 0(l/z)). (30) 



The second exponential in (pOD is dominating the asymptotic behaviour of ipi and ip2- 
Which one of these two will be the exponentially decreasing solution depends on the 
sign of the imaginary part of e'^*^^. As a consequence, the complex lower-half x-plane is 
divided into vertical stripes defined by 



Si = {x E C|Imx < 0,/7r < Rex < (/ + l)vr] , / G Z . 



(31) 



For X G 5*2^ we find that ipi is the normalizable solution. Whereas for x G 52;+! the 
normalizable solution is given by ip2- At the borders between neighbouring stripes, that 
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is, for Hex = In the first exponential in botli relations ( pOD guarantees square integrability 
of both solutions. 

The above discussion shows that the contour along which we have to define our quan- 
tum system has to be confined to one of these strips for asymptotically large x. In other 
words, the anti-Stokes lines for the present potential lie vertical in the lower-half plane 
and are equally spaced by a distance vr. In fact, this is what one expects in the limit 
cxD in (pSf). Hence, there are infinitely many possible choices in defining a proper 
quantum eigenvalue problem for the present potential. In what follows we will discuss 
some typical cases. 

Without loss of generality let us assume that the contour starts at Imx = — oo , say, 
in So- Hence, ipi is the proper solution. If we now choose our configuration space such 
that this contour also ends in the same sector we will not attain any restrictions on E. 
That is, we have again the unphysical situation of a quantum system with a spectrum 
unbound from below. Therefore, in order to have a proper quantum system we must 
demand that the contour ends in some other sector, say, Sm with m > 0. In this case we 
have to investigate the analytic continuation of ipi into this sector, which is given by [|13 



sin(z/7r) H^^^ {z e™^') = - sin [(m - 1)u7t] HI^'^ (z) - e'"^' sin(mz/7r) Z/^^) _ (32) 

where E = For m = 2/ > we have to avoid an admixture of H^'^\ which is the 

exponentially growing solution in 5*0 and S21. Hence we must demand that sin(2/7ri/) = 
and sin(7ri/) 7^ 0. In the particular case / = 1 this results in the condition u = n + ^, 
n = 0, 1, 2, 3, . . ., giving rise to the positive discrete spectrum 

En = l(n + l] , n e No . (33) 



2 V 2 



In the general case m = 21 we arrive at the conditions 2/z/ G N and u ^ N. For example, 

113 5 3 
4' 2' 4' 4' 2' 

1 

(2) 

Now we investigate the case m = 2Z + 1 for which Hu represents the normalizable 
solution in the corresponding sector. For the special case / = 0, that is, the contour ends 



for / = 2 we have u = leading to a "perforated" spectrum 



^" = ^(^^ + 1) ' riGNo\{2/-l,4/-l,...} . (34) 



in the neighbouring sector Si we have to invoke the connection formula ||T3 



Hi'\ze--^) = -e-^-^Hi'\z), (35) 

which does not yield any condition on u. In other words, the corresponding spectrum is 
unbound from below and thus unphysical. For / > 1 we again make use of (|3^) leading to 
2/z/ G N with z/ ^ N. These conditions are identically to those already found in the case 
m = 21. Hence, we will obtain the same spectrum (|34D . 

The remaining part of this section is devoted to a short semiclassical analysis of the 
problem (^). If we look at the potential at the left border of 5*0, that is Rex = 0, the 
potential has the form V{—iR) = |e^^ and thus exhibits a classical turning point for 
positive E. At the left border of 5*^, i.e. Rex = mn, the potential looks the same because 
of its periodicity. For a semiclassical analysis we need two complex classical turning points 
given, for example, by Xi = — iln/c and X2 = mn — iln /c, m > 1. They are the solutions 
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of fc^ — e^^^ = with fc^ = 2E. The semiclassical quantization condition is conventionally 
given by 

/ = ^ dx^/k^- e2i^ = n (^n + ^ , uENq. (36) 

In order to get a real value for the integral we have to integrate along the horizontal line 
between xi and X2 in the complex plane 0. Taking into account the periodicity of the 
integrant and changing the integration variable to 2; = e^^^/k'^ the integral can easily be 
calculated and yields / = mirk. Therefore the semiclassical approximation leads to 

1 / 1^^ 



^r = ^ , nGNo, (37) 



which, in fact, is very similar to the exact result (|3J) valid for both cases m = 21 and m = 
2Z + 1. However, we should note two essential failures of the semiclassical approximation. 
Firstly, it gives rise to a discrete spectrum even in the case m = 1 where the exact 
treatment has led to an unbound spectrum. Secondly, the semiclassical approach also 
fails to reproduce the perforation in the spectrum found for the cases m > 2, c.f. 
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Figure Captions 

Fig. 1: Real (a) and imaginary (b) part of potential (|12]) with (|16D for various values of e 
and fixed (3 = i. For purely imaginary (3 this potential is PT invariant. 

Fig. 2: Same as Figure 1 but for fixed e = —1/2 and various complex values of (3 with 
Re /3 = 0.5. Note that for these values of parameters the potential is not PT 
invariant. 

Fig. 3: Same as Figure 2 with Re /3 = 2. A singular behaviour emerges as P approaches 
the real axis. 
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